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ABSTRACT 


Context: Fast radio bursts are transient radio pulses of extragalactic origin. Their 
dispersion measure is indicative of the baryon content in the ionized intergalactic 
medium between the source and the observer. However, inference using unlocalized 
fast radio bursts is degenerate to the distribution of redshifts of host galaxies. 

Method: We perform a joint inference of the intergalactic baryon content and the 
fast radio burst redshift distribution with the use of Bayesian statistics by comparing 
the likelihood of different models to reproduce the observed statistics in order to infer 
the most likely models. In addition to two models of the intergalactic medium, we 
consider contributions from the local environment of the source, assumed to be a 
magnetar, as well as a representative ensemble of host and intervening galaxies. 

Results: Assuming that the missing baryons reside in the ionized intergalactic 
medium, our results suggest that the redshift distribution of observed fast radio bursts 
peaks at z < 0.6. However, conclusions from different instruments regarding the inter- 
galactic baryon content diverge and thus require additional changes to the observed 
distribution of host redshifts, beyond those caused by telescope selection effects. 


Key words: cosmology: observations — cosmology: large-scale structure of universe — 
galaxies: intergalactic medium — galaxies: distances and redshifts — radio continuum: 


transients — transients: fast radio bursts 


1 INTRODUCTION 


Observations of the cosmic microwave background (CMB) 
show that ~ 5% of energy density in the Universe exists 
as ordinary matter (Planck Collaboration et al. 2014). 
However, in the z < 2 range, observations of stars and 
gas in galaxies, of the hot intracluster medium and of the 
Lya-forest can only account for about half of that amount 
(e. g. Nicastro et al. 2008). The ”missing” baryons are 
believed to be hidden in the warm-hot intergalactic medium 
(WHIM) (Cen & Ostriker 1999; Dave et al. 2001), with 
temperatures of 10° — 10’ K and low baryon densities of 
107 — 107" cm~’. This medium is difficult to observe 
directly as it only couples with radiation through electronic 
transitions. However, correlation of the distribution of 
galaxies with CMB distortions due to Sunyaev-Zel’dovich 
effect has revealed the likely presence of filaments of 
warm-hot gas between galaxies, possibly consistent with 
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the missing baryons (e. g. Tanimura et al. 2018; de Graaff 
et al. 2019). 


Several studies have suggested that fast radio bursts 
(FRBs) can be used to detect the missing baryons. Their 
signals propagate across cosmological distances and get dis- 
persed by diffuse ionized gas, quantified by the dispersion 
measure (DM), defined as the column density of free elec- 
trons 


DM = | nedi. (1) 


Wei et al. (2019) proposed to compare the redshift evolution 
of the DMs using ~ 3000 FRBs with known redshifts to 
that of the Hubble parameter to measure the IGM baryon 
content fram, i. e. the amount of all baryons that reside 
in the ionized IGM. Deng & Zhang (2014) argued that 
the redshift observed for y-ray bursts associated to FRBs 
can be used to infer fram. McQuinn (2013) modelled 
the distribution of DM expected for FRBs from different 
redshift and found that ~ 100 DM of FRBs from redshift 
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z > 0.5, localized with sub-arcminute precision around 
identified galaxies, can be used to infer the baryon profile 
of galaxies in order to constrain the localization of cosmic 
baryons. Muñoz & Loeb (2018) proposed to cross-correlate 
the DM of arcminute-localized FRBs with the thermal 
Sunyaev-Zeldovich effect in the CMB. The latter depends 
on the temperature of the WHIM and the amount of 
baryons localized in the WHIM, which can be constrained 
using ~ 1000 FRBs with arcminute localization. Qiang & 
Wei (2020) suggest to use model-independent Gaussian 
processes to investigate the evolution of firem with DM of 
FRBs with identified redshift. 


All of the aforementioned papers require a large 
number of well-localized FRBs or cross-correlation with 
associated observations. However, Yang et al. (2020) 
argue that the majority of FRBs are not supposed to 
have strong associated persistent sources. Furthermore, 
the exact localization of sources of short-duration signals 
without known redshift is not trivial to constrain, and 
neither is the identification of the host galaxies of FRBs 
(Eftekhari & Berger 2017; Mahony et al. 2018; Marcote 
& Paragi 2019; Prochaska et al. 2019a). The host can be 
identified using interferometry or by observing persistent 
counter parts to the FRB. So far, FRBs have been located 
to a multitude of galaxy types, from star-forming dwarf 
galaxies (Tendulkar et al. 2017) to very massive galaxies 
with old stellar population (Ravi et al. 2019; Bannister 
et al. 2019). Unluckily, the small sample of localized FRBs 
cannot provide enough information to arrive at reasonable 
conclusions regarding cosmological questions. However, in 
this work, we show how it is possible to use DM of a large 
sample of unlocalized FRBs, in order to infer their host 
redshift distribution in a statistical way, which also allows 
us to constrain the amount of ionized baryons located in 
the cosmic web. 


Because of the unknown distance to unlocalized FRBs, 
the assumed distribution of source redshifts can help inter- 
pret the distribution of observed DM. Several papers have 
tried to infer the intrinsic redshift distribution of FRBs, ei- 
ther by modelling the distribution of DM and other FRB 
properties with analytical or Monte-Carlo methods (Bera 
et al. 2016; Caleb et al. 2016; Gardenier et al. 2019), or by 
performing a luminosity-volume test (Locatelli et al. 2018). 
So far there are inconsistencies in implications of data from 
different instruments, e. g. ASKAP requiring faster change 
in FRB density with redshift than Parkes. 


Macquart & Ekers (2018) use the flux density and 
fluence of FRBs to infer the luminosity distance and 
evolutionary history, as well as the redshift distribution, 
and show how the history of the ionized IGM affects the 
distribution of observed DM, thus influences our inference 
of the host redshift distribution. 


Here we propose to use unlocalized FRBs to perform 
a joint analysis of the FRB redshift distribution and the 
IGM baryon content, firem, by comparing the expected dis- 
tribution of DM to the available observed values reported by 
Parkes, CHIME and ASKAP observatories (Staveley-Smith 
et al. 1996; Macquart et al. 2010; The CHIME/FRB Collab- 


oration et al. 2018), using the PREFRBLE' software (Hack- 
stein 2020), presented in Hackstein et al. (2020). In Sec. 2 
we explain how we model different values of fiam. The re- 
sulting predictions are presented and compared in Sec. 3. 
We discuss our results in Sec. 4 and conclude in Sec. 5. 


2 METHOD 


Investigation of the fraction of baryons in the ionized IGM, 
fiam = Niam /Nbvaryons, With DM requires comparison of the 
observed DM with expectations of DM(z|fiam), according 
to source redshift z, (Keane et al. 2016) 


CPcrits 4b Qe (1 a z ) 
(DM(2)) = Pee f free) Gary” ae (2) 
However, FRBs do not come with a direct measure of red- 
shift and for the majority of FRBs, that could not be local- 
ized by other means, the DM is the best indicator for the 
source distance (Dolag et al. 2015; Niino 2018; Luo et al. 
2018; Walker et al. 2018; Pol et al. 2019). Still, by assum- 
ing different plausible redshift distributions of FRBs, x(z), 
we can estimate the distribution of DM to be observed by 
instruments, which is determined by fiam. 

By comparing the expected distribution to observations, 
we can quantify the likelihood of different combinations 
of m(z) and fiam. Consequently, this allows us to put 
constraints on the WHIM density without the need to 
localize FRBs. Furthermore, by a detailed investigation of 
the redshift-evolution of fram in cosmological simulations, 
the inference presented here can be used to investigate the 
helium-reionization history (e. g. Linder 2020; Dai & Xia 
2020). 


In Hackstein et al. (2020), we obtained likelihood esti- 
mates, L(DM|z), to observe extragalactic DM from source 
at redshift z. These expectations consider contributions from 
all regions along the LoS in our benchmark scenario, consid- 
ering magnetars as the source of FRBs (Piro & Gaensler 
2018) as well as a realistic ensemble of different host and 
intervening galaxies (Lacey et al. 2016; Rodrigues et al. 
2018). Predictions for the IGM were obtained using a con- 
strained cosmological simulation of the local Universe, pro- 
duced using the cosmological magneto-hydrodynamic code 
ENZO (Bryan et al. 2014) together with initial conditions 
obtained following Sorce et al. (2016) and cosmological pa- 
rameters of PLANCK (Ade et al. 2016). The constrained 
volume of 250 (Mpc/h)? was embedded in a 500 (Mpc/h)? 
grid of 1024° cells. The use of five adaptive mesh refinement 
levels allowed for a maximum resolution of 30 kpc in the 
most-dense environments. The simulation starts at redshift 
z = 60, where all baryons are in the IGM. However, the lim- 
ited resolution does not allow us to properly resolve galaxy 
formation and the condensation of cold gas out of the IGM, 
thus implies fiam = 1, always. Further information on this 
model can be found in Hackstein et al. (2018) and Hack- 
stein et al. (2019). A reduced version of this model can be 
found on crpropa.desy.de under “additional resources”, to- 
gether with the other models presented in Hackstein et al. 
(2018). 


1 sithub.com/FRBs/PreFRBLE 


MNRAS 000, 1-9 (0000) 


The contribution of the circumburst environment of 
magnetars is modelled using the analytical predictions of 
the wind model in Piro & Gaensler (2018) in a Monte-Carlo 
simulation. More details on this simulation and the assumed 
priors can be found in Hackstein et al. (2019). The contri- 
bution of host and intervening galaxies is estimated by com- 
puting LoS integrals in a realistic ensemble of galaxy models 
that recreate a wide variety of properties of the population of 
galaxies in the Universe (Lacey et al. 2016; Rodrigues et al. 
2018). More details on this Monte-Carlo simulation and the 
resulting likelihoods can be found in Hackstein et al. (2020). 

We further use FRBPOPPY” (Gardenier et al. 2019) 
to model the distribution of host redshift of observed 
FRBs, m(z). We assume three physically motivated intrin- 
sic redshift distributions, i. e. following the comoving vol- 
ume (coV), stellar mass density (SMD) or star formation 
rate (SFR), together with the selection effects of ASKAP, 
CHIME and Parkes instruments. We use the parameters of 
their complex model, which best reproduces observations. 

Together, these can be used to quantify the distribution 
of expected DMs, 


L(DM) = [rome x m(z) dz. (3) 


For more details on the derivation of likelihood L(DM|z) 
and prior 7(z), see Hackstein et al. (2020). 


Numerical simulations suggest that up to 28 + 11 per 
cent of baryons are located in neutral hydrogen clouds in the 
IGM, observable in the photoionized Lya forest, as well as 
25 + 8 per cent in the WHIM, leaving a baryon shortfall of 
29 + 12 per cent (Shull et al. 2012). This leaves a plausible 
range of 0.3 S fram £ 0.9 (cf. Li et al. 2019). 

As explained above, in the constrained simulation friem 
is assumed to be constant with redshift, thus to be a global 
factor to (DM(z)) (Eq. 2). This allows to obtain estimates 
for different values of frcm in post-processing, allowing us 
to test how many missing baryons can be located in the 
WHIM (Deng & Zhang 2014; Keane et al. 2016). To this 
end, we modify the likelihood of intergalactic contribution, 
L(DMiem|z), obtained from the constrained simulation of 
the IGM, such that the result returns the required average 
value according to Eq. 2, while shape and norm are con- 
served: 


L(DMiem|fiem) = fiam X L(fiam x DMiem). (4) 


The fraction of baryons in the IGM, frcm, is believed 
to grow with redshift (McQuinn 2013; Prochaska & Zheng 
2019). Simulations suggest that fram % 0.9 prior to z = 1.5 
(Meiksin 2009). To account for the redshift evolution of 
fiam, we consider a model for the IGM following Pshirkov 
et al. (2016) (called Pshirkov16 hereafter). We perform 
a Monte-Carlo simulation in which we divide the LoS in 
segments of Jeans-length size \;(z) = 2.3 (1+ z)~'° Mpc. 
Assuming that the electron density ne is well described by 
the distribution of neutral hydrogen deduced from the ob- 
served Lyman-a forest, we pick random values for comoving 
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ne drawn from a log-normal distribution of overdensity 
according to Eq. (2) in Pshirkov et al. (2016). Finally, ne is 
multiplied by the value of fiam at corresponding redshift. 
We assume fram to increase linearly with redshift from 
ficam = fiam,o at z = 0 to friem = 0.9 at z = 1.5, while it 
stays at this constant value at higher redshifts. 


In order to quantify the likelihood of a scenario of friem 
and FRB redshift distribution 1(z) to account for the obser- 
vation of an individual FRB with DMo»s = DMea+DMmw, 
we take the value of L(DMrzc| ficam) for the benchmark sce- 
nario, using the estimates for DMmw according to NE2001 
model (Cordes & Lazio 2002) as listed in FRBcat (Petroff 
et al. 2016). These estimated values of DMrc were shown 
to be correct to & 30 pe cm~* (Manchester et al. 2005). 

In order to compare different choices of fiam and host 
redshift distribution, we compute the Bayes factor 


L(DM i 
B(DMrc|ficm,i, fIGM,ref) = Ai L 


by renormalizing all likelihoods of individual FRBs to the 
same reference choice of fiGm,ref. B quantifies how much the 
observation of DMgeq changes our corroboration in favour 
of fiam,; compared to figm,rer- By renormalizing all B to 
the best-fit scenario Bmax, the presented values quantify how 
much the observations change our corroboration in disfavour 
of the corresponding scenario. 

For FRBs with localized hosts, the information on 
fiam and z(z) can be distinguished, thus allowing for 
separate likelihoods. For the events listed with host redshift 
zFRB in FRBcat, instead of the procedure explained above, 
we obtain the likelihood of firem by comparing values of 
L(DMge|zrrp, ficm) and compute the Bayes factor B 
by renormalizing the likelihood to a particular choice of 
fiam,. Likewise, we compare values of 7(zpre) for their 
implication on the redshift distribution. 


(5) 


B obtained for individual events and different instru- 
ments can simply be multiplied in order to interpret the 
combined results. Note that 6B quantifies how much the cur- 
rent observations change our corroboration between the con- 
sidered scenarios. However, according to Bayes theorem 


P(ficm|DMrc) x L(DMge|fiem)7(frem), (6) 


the ratio of likelihoods L, 6, has to be multiplied by the 
ratio of priors m in order to arrive at the ratio of posteri- 
ors P. n( fiam) quantifies how likely different values of frem 
are prior to investigation of DMra, while P expresses our 
conclusion on the likelihood of values of fiam after consid- 
eration of the data. However, we consider no preference for 
any value of fiam, 7(ftam) = const., thus B is identical to 
the ratio of P. 


3 RESULTS 
3.1 Expected distribution vs. observation 


In Fig. 1 we show the results for all combinations of instru- 
ments and redshift distributions for two extreme constant 
values of frem = 0.3 and frem = 0.9, using the constrained 
simulation for the IGM. These can be compared directly 
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constrained simulation, friem = const. 
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Figure 1. Distribution of extragalactic DMge = DMops — DMmw as observed (dash-dot) by ASKAP (left), CHIME (center) and Parkes 
(right) minus estimate of MW contribution based on NE2001 model, as listed in FRBcat. The solid lines show the expected distribution 
to be observed by the corresponding instrument according to Eq. 3, assuming FRBs from magnetars in our benchmark scenario using the 
constrained simulation to model the IGM and extreme constant values for the baryon content, fram = 0.3 (top) or fram = 0.9 (bottom), 
as well as FRB redshift distribution to follow star formation rate (SFR, blue), comoving volume (coV, orange) or stellar mass density 
(SMD, green). The error bars show the shot noise of the observed data. For the expected distribution, the barely visible error bars are 
the shot noise of the Monte-Carlo samples used to obtain the likelihood function. The product of Bayes factors (Eq. 5) for individual 
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FRBs, shown in Fig. 3, quantifies how well different scenarios recreate the observed distribution. 
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Figure 2. Same as Fig. 1, using the Pshirkov16 model with redshift dependent fragm for the IGM and extreme values for its present value 
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ftam,o- Models for the redshift distribution, circumburst environment, host and intervening galaxies remain the same. 
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to observations by the corresponding instrument, which are 
indicated by the dash-dotted line. 

The majority of FRBs observed by ASKAP have 
DMgc & 500 pe cm™?. The expected peak values of DMga 
for fram between 0.3 and 0.9 range from 100 pc cm~® to 
440 pc cm~*. ASKAP results are thus in favour of high val- 
ues of fiam % 0.9. However, for the values of DM observed 
with ASKAP that are generally higher than expectations, 
other foregrounds such as the local contributions need to be 
considered more carefully in order to arrive at reasonable 
conclusions. The assumed intrinsic redshift distributions do 
not differ much in their predictions at low z (cf. Hackstein 
et al. 2020), thus rendering ASKAP results incapable of dis- 
tinguishing FRB populations with a small sample of FRBs. 

While the distribution of DM observed by CHIME 
peaks at roughly the same value as for ASKAP DMge ~% 
500 pc cm~?, the peak is less pronounced and the upper 
1o-deviation is more than 1.5 times higher than for ASKAP, 
reaching values beyond 1000 pc cm~?. However, for CHIME, 
the expected peak values of DMgg range from 300 pc cm~? 
to 1500 pc cm~® for fram = 0.3 and fram = 0.9, respec- 
tively, far above the observed value, hence expecting many 
more high values of DM. Thus, the CHIME results favour 
low friem Z 0.3, depending on the assumed distribution of 
the host redshifts. Also, CHIME observes more FRBs with 
50 < DMgc / pc cm~? < 200 than ASKAP, which is not 
expected according to results in Fig. 1. This hints on an 
imprecise assumption on the redshift distribution of FRB 
sources visible to CHIME, that observes at different wave- 
lengths than the other two instruments. 

For Parkes, the expected peak value of DMre ranges 
from 300 pc cm~? to 1800 pe cm™? between friem = 0.3 
and frem = 0.9, similar as expected for CHIME. However, 
the observed sample peaks at DMga œ 1000 pc cm~? - 
about twice the value observed with CHIME - which is 
well within the expected range. Hence, different assumed 
redshift distributions will favour other values of friem, e. g. 
a population concentrated at low redshifts, such as SMD, 
will need a high value of firem < 0.9 to explain the observed 
DMs, while populations with greater redshifts, such as 
SFR, will require lower values of frem. The Parkes sample 
is dominated by high values of DMs, with ~ 30 per cent of 
FRBs observed with DM > 1000 pc cm~?. The majority 
of these is thus likely located at z 2 1, such that local 
contributions to the DMs are less significant. This makes 
the Parkes results overall the most reliable for investigation 
of friem and other cosmological properties. 


In Fig. 2 we present the results for all combinations 
of instruments and redshift distributions for two extreme 
values of fiam,o = 0.3 and fiam,o = 0.9 at z = 0, by mod- 
elling the IGM with the Pshirkov16 model, including the 
linearly increasing friem. For the case that fiam,o = 0.9, 
Jiem = const. at all redshifts, similar to the constrained 
simulation in the lower row of Fig. 1. The fact that predic- 
tions are hardly distinguishable verifies the similarity of both 
models ? However, for lower values of fremo there is stronger 
contribution for DM at large redshift, thus larger values of 
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DM are predicted in this case. This reduces the variance of 
predictions with different ftam,o, when compared to varying 
constant fiam. 


3.2 Inference of frcm and 7(z) 
8.2.1 fiam from localized FRBs 


Recently, Li et al. (2020) estimated fiam = 0.847938 us- 
ing DM of the five localized FRBs together with the cor- 
responding luminosity distance of Ia supernovae. Shortly 
after that, Bhandari et al. (2020) report the localization 
of two additional FRBs. Using their reported values of 
DMec = DMops — DMmw for all six FRBs localized at de- 
cent redshift, z > 0.1, we estimate figm = 0.8470 35 within 
30 using our more sophisticated approach with our bench- 
mark scenario and the Pshirkov16 model for the IGM. The 
lower limit suggests that friem Z 0.5. The current sample of 
localized FRBs is thus not yet sufficient to solve the missing 
baryons problem unambiguously (cf. Macquart et al. 2020). 

Though our inference depends on the chosen cosmology 
(Ade et al. 2016), as the models for IGM as well as host 
and intervening galaxies in our benchmark scenario are ob- 
tained from a specific cosmological simulation, we derive the 
same estimated value, with a more conservative lower limit. 
Still, the conclusion that this confirms the location of miss- 
ing baryons in the WHIM should be taken with caution, at 
least because frem is considered to be constant with red- 
shift. Note that the likelihood function used in this estimate 
for fiom, i. e. [] L(friem|DM;, zi), can be used as prior to 


infer the host redshift distribution using unlocalized FRBs, 
replacing assumptions on frcm. 


3.2.2 fiam from individual instruments 


In Fig. 3 we show the Bayes factors B resulting from the joint 
analysis of the FRB redshift distribution and of the IGM 
baryon content fram, assumed to be constant in our con- 
strained simulation. For the 28 FRBs observed by ASKAP, 
high values of 0.8 < fiam are strongly favoured, regardless 
of the assumed redshift distribution. Hence, ASKAP val- 
ues seem to confirm the location of missing baryons in the 
ionized IGM, most probably the WHIM. Due to the gen- 
erally low distance of ASKAP events, B shows the lowest 
variability between the assumed host redshift distributions. 
However, estimates using ASKAP results are more vulnera- 
ble to the unknown contribution of the host galaxy and local 
environment of the source and should be taken with great 
caution. This is further stressed by the fact that ASKAP 
likely prefers too high friem Z 0.9. 

The 30 FRBs observed by CHIME, of which 18 have 
been found to repeat (CHIME/FRB Collaboration et al. 
2019; Andersen et al. 2019; Fonseca et al. 2020), clearly 
favour low values of frem < 0.5 as the observed distribution 
of DM concentrates on lower values than expected. Since 
SFR redshift distribution peaks at a higher redshift, z ~ 1, 
it is disfavoured by CHIME results, except for the case of 
fiam ~ 0.3. Instead, the SMD distribution, which peaks at 
lowest redshift, z ~ 0.5 is preferred for all probed values of 
ficm > 0.4 and favours friem = 0.4. Only for firem < 0.4, 
coV and SFR show higher values of B. However, the param- 
eters for CHIME selection effects available in FRBPOPPY are 
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constrained simulation, friem = const. 
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Figure 3. Bayes factor B of FRBs found in FRBcat comparing scenarios for different baryon content in IGM frcm and redshift distributions 
of FRBs, following either star formation rate (SFR), comoving volume (coV) or stellar mass density (SMD), assuming that FRBs are 
produced by magnetars in our benchmark scenario. To model the IGM we use the constrained simulation described in detail in Hackstein 
et al. (2020). The title indicates the instrument that observed the considered number of N FRBs. The estimated value and deviation 
are obtained via Jackknife resampling. All B are renormalized to the best-fit Bmax to indicate how much less likely a given scenario 


reproduces the observed data (cf. Fig. 1). 
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Figure 4. Same as Fig. 3, using the Pshirkovl6 model for the IGM contribution. Models for the redshift distribution, circumburst 


environment, host and intervening galaxies remain the same. 


an early version (see Gardenier et al. 2019), hence results for 
this instrument should be considered with caution. Still, if 
the current parameters are proven to be at least roughly cor- 
rect, the present results either indicate a wrong assumption 
entering the expected distribution of FRB redshifts, which 
is likely as CHIME observes a different frequency band than 
ASKAP and Parkes, or suggest a very low baryon content 
in the IGM, which we deem as unplausible, as it would con- 
tradict the result of any cosmological simulation. However, 
CHIME results so far seem to disfavour the "WHIM solu- 
tion” to the missing baryons problem. 


The 27 FRBs observed by Parkes show a high variation 
among the assumed host redshift distributions. The SFR 
redshift distribution peaks at highest redshift z ~ 1.1 and 
favours 0.4 < fiam < 0.7, while it cannot be brought into 
agreement with the WHIM solution, as B(fiam > 0.7) < 
107? indicates that these values can be ruled out with more 
than 99 per cent certainty. The SMD distribution, peaking 
at z © 0.6 prefers a higher fram > 0.7 value. The coV dis- 
tribution peaks at z ~ 0.7 and favours friem œ% 0.6, however, 
with relativley low significance. 


In Fig. 4 we show the Bayes factors B resulting from 
the joint analysis of the FRB redshift distribution and of the 
IGM baryon content fiam,o at z = 0, assumed to increase 
linearly until firem = 0.9 at z = 1.5 in the Pshirkov16 
model. Since predictions show less variance in fiam,o, the 
graphs are shallower than in Fig. 3, while the general 
trends remain the same. ASKAP results are in favour of the 
WHIM solution to the missing baryon problem. However, 


CHIME results strongly disfavour high densities of baryons 
in the local IGM, for all assumed redshift distributions. In 
particular, there is no value of ftam,o for which B > 107? 
with both, ASKAP and CHIME. Such results can thus 
not be interpreted consistently, without considering further 
features that differ between the telescopes, e. g. systematic 
differences in the local contribution or the population of 
potentially visible FRBs. However, investigation of possible 
solutions is beyond the scope of this paper. 


Not only does the inference of the different data sets 
differ dramatically, but also the assumed redshift distribu- 
tion of FRBs strongly affects our interpretation of the data. 
Thus, the degenerate problem of inferring friem with unlo- 
calized FRBs is unlikely to be solved simply by increasing 
the sample size. This stresses the importance of careful in- 
vestigation of all assumptions entering the interpretation of 
measures of FRBs in order to arrive at reasonable conclu- 
sions. As for now, we cannot derive strong or unambiguous 
conclusions on frem and on the intrinsic redshift distribu- 
tion from unlocalized FRBs. Only for the CHIME results 
we can conclude that the underlying redshift distribution 
has to peak at z < 0.6, in order to explain the high num- 
ber of sources with low DMs. This conclusion is independent 
of the details of the assumed intrinsic redshift distributions 
and selection effects. 
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4 DISCUSSION 


Our work highlights the fact that the distributions of FRBs 
obtained by different radio telescopes lead to very different 
inferences about the cosmological distribution of such 
events. The disagreement between ASKAP and CHIME is 
to be explained by a strong local contribution of DM to 
the bright low-distance FRBs observed by ASKAP. These 
would account for systematically increased values of the 
observed DM. By not accounting for such a strong local 
contribution in probability estimates, the increase in DM 
would be misinterpreted as IGM contribution, calling for a 
high fram. If the high DMs observed with ASKAP can be 
attributed to a stronger local contribution, lower values of 
fiam would be favoured by FRB observations. This would 
reduce the fraction of missing baryons found in the WHIM. 


However, the fact that CHIME sees a higher amount 
of FRBs with 50 < DM / pecm~* < 200 than ASKAP 
cannot be explained by low baryon content in the IGM or the 
general FRB population to be concentrated at low redshifts, 
which would affect CHIME and ASKAP in the same way. If 
this feature persists in larger samples, it provides a strong 
hint towards CHIME observing a different subset of FRB 
sources than ASKAP. The more than 10 times higher gain 
of CHIME allows it to potentially measure a different set of 
fainter sources (Amiri et al. 2018), e. g. magnetars of older 
age (cf. Metzger et al. 2017; Marcote et al. 2020). Since 
these fainter sources are only visible at low redshift, this 
would increase the likelihood for CHIME to observe FRBs 
with low DM. Furthermore, the local environment of older 
magnetars contributes much less to observed DM (Piro & 
Gaensler 2018). Thus, old magnetars could account for the 
increased number of FRBs observed with low DM compared 
to our expectations that consider the same population of 
magnetars to be observed by CHIME and ASKAP. Such an 
increase in the likelihood to observe FRBs from low redshift 
would also make the CHIME results favour higher values of 
fiam. However, since CHIME is also able to observe FRBs 
visible with Parkes (e. g. FRB121102, see Josephy et al. 
2019), we require CHIME to detect more events with DM > 
800 pc cm~? and any above > 2000 pc cm~® in the future. 
If there is an extended set of sources available, this should 
reflect in a higher rate of FRBs observed by CHIME than 
extrapolated from other surveys (e. g. Chawla et al. 2017). 


Note that 18 of the 30 FRBs observed by CHIME we 
considered here are repeaters (CHIME/FRB Collaboration 
et al. 2019; Andersen et al. 2019; Fonseca et al. 2020). The 
high number of repeaters already observed by CHIME would 
suggest an increased chance to observe fainter bursts from 
repeating sources, e. g. from old magnetars. For these young 
neutron stars, this implies a more advanced dissipation of 
supernova remnants, thus a lower local contribution to DM. 
If ASKAP can only observe magnetars in an early stage, 
say < 1000 yr, this could account for the high local contri- 
bution to DM, required above, but would also lower the in- 
ferred value of fiam. Older magnetars potentially observed 
by CHIME would also cause a lower contribution to RM, 
which could be identified by statistical investigation of the 
distribution of RM observed by the different instruments. 
However, the current sample of reported RM is too small 
to allow for any firm conclusions. It might also be possible 
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that separate frequency bands are dominated by completely 
different sources of FRBs. 

Still, even with negligible local contributions, the 
CHIME results suggest a low IGM baryon content friem S 
0.5. If this cannot be explained by false assumptions about 
the FRB redshift distribution and selection effects, the 
CHIME results would constitute strong evidence against 
missing baryons in the WHIM or, in fact, anywhere in the 


ionized IGM. 


Caveats Our conclusions are based on a number of 
assumptions regarding the IGM baryon content, the FRB 
redshift distribution as well as local contributions and 
foregrounds. Furthermore, in our benchmark scenario we 
assumed all FRBs to be caused by magnetars in axisym- 
metric galaxies. Alternative progenitor models beyond our 
assumptions might lead to different conclusions. For exam- 
ple, magnetars that only produce FRBs at an early stage, 
say less than 10° yr after their formation, would produce 
much stronger dispersion than older ones (> 10° yr), which 
we allowed in our benchmark scenario. This would make 
lower values of fiam more plausible, especially for the 
low-distance sample observed by ASKAP. In future work, 
we will consider a wider range of assumptions in order to 
provide more significant results. 


Note that FRBPOPPY uses estimates, e.g. of DM(z), 
in order to decide how many FRBs will be observed at 
a given redshift. Theses estimates have been produced 
using slightly different assumptions on the contributing 
regions. However, the DM is dominated by the IGM and 
the analytical description used in FRBPOPPY provides a 
good match to our estimates. Hence, we argue that this 
does not alter the general conclusions of this work. 


Finally, in our benchmark scenario, we do not account 
for clustering or galaxy haloes (e. g. Prochaska et al. 2019b; 
Connor et al. 2020). These are especially important for the 
interpretation of FRBs located at low redshifts. However, 
FRBs from higher redshift, for which these contributions are 
less important, are much more informative regarding 7(z) 
as well as frem. Our results are therefore reasonably robust 
against the additional presence of haloes along the line of 
sight. 


5 CONCLUSIONS 


Fast Radio Bursts are an important probe of Baryons in 
the Universe, as their dispersion measure is sensitive to all 
phases of the intergalactic medium, including the WHIM. 
However, inference of the location of missing baryons with 
unlocalized FRBs is degenerate to the distribution of their 
host redshifts. In this work we use the estimated extragalac- 
tic component of the dispersion measure to jointly infer the 
redshift distribution of unlocalized FRBs and of the IGM 
baryon content with a Bayesian statistical approach. We use 
an analytical model for the local environment of the source, 
assumed to be a magnetar, semi-analytic models for the en- 
semble of host and intervening galaxies, as well as models for 
the IGM from cosmological simulations. Comparing expec- 
tations to observations of unlocalized FRBs, we investigate 
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the implications of events observed by ASKAP, CHIME and 
Parkes, listed in the FRBcat, using the open-source python 
software package PREFRBLE* (Hackstein 2020). Our main 
conclusions can be so summarised: 


e From six localized FRBs beyond redshift z > 0.1, as 
listed in Bhandari et al. (2020), we infer the intergalactic 
baryon content with 30 limits friem = 0.844911, FRBs thus 
potentially confirm the location of missing baryons in the 
WHIM. This value agrees well with results of Li et al. (2020). 
However, the error margin does not yet allow to confirm the 
WHIM solution beyond doubt. 

e In order to unambiguously infer frem from unlocalized 
FRBs, the discrepancy between different instruments need 
to be resolved, e. g. by identifying different subsets of the 
FRB population. A mere increase in the number of events 
will likely not solve this discrepancy. 

e In case that the WHIM solution to the missing baryon 
problem is correct, the Parkes sample is in conflict with a 
FRB population following the star formation. 

e Limited to the CHIME sample, we find that the distri- 
bution which tracks the star formation rate is the least likely 
to explain observations. A better agreement is achieved for 
a distribution that peaks at rather low redshift z < 0.6, e. 
g. following the stellar mass density. 


Though the small sample of localized FRBs suggest the 
missing baryons are in the ionized IGM, the error margin 
does not yet allow for unambiguous confirmation of this hy- 
pothesis. 
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